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ABSTRACT 

We present new Ryle Telescope (RT) observations of the Sunyaev-Zel'dovich (SZ) 
decrement from the cluster Abell 773. The field contains a number of faint radio 
sources that required careful subtraction. We use ASCA observations to measure 
the gas temperature and a ROSAT HRI image to model the gas density distri- 
bution. Normalizing the gas distribution to fit the RT visibilities returns a value 
of Hq of 77l|^5 kms~^ Mpc~^ (l-tr errors) for an Einstein-de-Sitter universe, or 
85ti?kms-iMpc-^ for a flat model with Qa ^ 0.7. The errors quoted include es- 
timates of the effects of the principal errors: noise in the SZ measurement, gas tem- 
perature uncertainty, and line-of-sight depth uncertainty. 
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1 INTRODUCTION 

We have previously reported the detection of a Sunyaev- 
Zel'dovich (SZ) decrement (Sunyaev & Zel'dovich 1972) to- 
wards the 2 = 0.217 cluster Abell 773 using the Ryle Tele- 
scope (RT) (Grainge et al 1993). (The SZ effect in this clus- 
ter has also been mapped by the millimeter array of the 
Owens Valley Radio Observatory (Carlstrom, Joy & Grego 
1996).) The RT observations of AbeU 773 form part of a con- 
tinuing programme to observe an X-ray luminosity-limited 
sample of rich, intermediate-redshift clusters in order to 
measure Ho by combining SZ and X-ray observations (Jones 
et al.2001). Such programmes (e.g. Reese et al., 2002; Mason, 
Myers and Readhead, 2001; see also Birkinshaw, 1999 for a 
review) are direct measurements of Ho free from distance- 
ladder arguments. 

In Grainge et al. 1993 we did not calculate an estimate 
of Ho because no suitable X-ray image of A773 and no esti- 



mate of its gas temperature existed. A ROSAT HRI image 
and ASCA spectroscopic data have since become available, 
and we have also made additional RT observations. These 
now enable us to make an estimate of the Hubble constant 
from this cluster, which, when combined with other clusters 
from the sample, will give an estimate of Ho unbiased by 
the individual shapes and orientations of the clusters. 



2 RYLE TELESCOPE OBSERVATIONS AND 
SOURCE SUBTRACTION 

The RT (Jones 1991) is an east-west synthesis telescope of 
13-m antennas with a bandwidth of 350 MHz and an aver- 
age system temperature for these observations of 65 K at an 
observing frequency of 15.4 GHz. We used five antennas in 
a compact configuration, giving two baselines of 18 m, three 
of 36 m, and five more out to 108 m. The short baselines 
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alone are sensitive to the SZ signal; the longer ones arc used 
to recognize and subtract the radio sources in the field that 
would otherwise mask the SZ decrement. We have made a 
total of 30 12-h observations of A773, each with the point- 
ing centre RA 09*" 17™ 51". 91, Dec. +51° 43' 32" (J2000). 
Phase calibration using 0859+470 and flux calibration using 
3C 48 and 3C 286 were carried at as described in Grainge 
et al (1993). Similarly, we used the POSTMORTEM package 
(Titterington 1991) to flag the data for interference and an- 
tenna pointing errors, and to weight them in accord with 
the continuously monitored system temperature of each an- 
tenna. As a standard check, we used the AiPS package to 
make a map of each 12-h run and then combined the data. 

We removed radio sources from the data by a simulta- 
neous maximum-likelihood fit to several point sources and 
the SZ effect using a technique described by Grainger et al 
(2002). We use a model for the SZ signal as a function of 
baseline that is based on the /9-model fit to the X-ray im- 
age described below (Section 3). Wo simultaneously fit flux 
densities for trial sources whose initial positions arc deter- 
minied both from a map made from just the long-baseline 
data (> 2 kA), and from a VLA 1.4-GHz image of the cluster 
field (Figure 1). This allows us to fit the optimum fiux den- 
sities of sources whose existence we know of from the VLA 
image but which would not give a significant detection from 
the RT data alone. The postitions and fitted flux densities 
are given in Table 1. The image made from the long (> 2kA) 
source-subtracted baselines is consistent with noise (Figure 
2). 

To image the decrement, we removed the sources in Ta- 
ble 1 from all the visibilities and made a short-baseline map 
from baselines shorter than 1 kA, and CLEANed this. The 

resulting image is shown in Figure 3. The docromont is of 
— 527/i Jybeam"^ with a noise (1-cr) of 60/iJybeam~^; the 
beam is 152 x 119 arcsec FWHM. Also shown is the X-ray 
image of the cluster; it can be seen that the alignment with 
the X-ray image is very good. The extension of the SZ image 
to the north-east is of marginal signiflcance. The magnitude 
of the decrement is consistent with that of —590 ± 116 /xJy, 
in the same beam, reported in Grainge et al (1993). 

An alternative way of looking at the data is shown in 
Figure 4, which shows the real part of the source-subtracted 
visibilities binned radially, along with the best-fitting model 
based on the X-ray data. These data have the advantage, un- 
like the image pixels, of having independent gaussian noise 
on each point; it is these that are used in the fitting for Hq. 



3 X-RAY OBSERVATIONS AND FITTING 

Wo measure the gas temperature from ASCA observations 
on 1994 April 29 of 46240 s (GIS) and 39904 s (SIS), using 
standard XSPEC tools. Times of high background flux were 
excluded and both GIS and SIS data were used. We took the 
Galactic absorbing column density predicted by Dickey and 
Lockman (1990) in the direction of A773 of 1.3 x 10^* H 
atoms m~^. Using a Raymond-Smith model, we find a tem- 
perature of 8.7±0.7 keV (90%-confidence error bounds) and 
a metallicity of 0.25 solar. The 2 10 keV flux from A773 is 
6.7±1.0 X 10^^'^ W m^^. Our temperature estimate is con- 
sistent with that of Allen and Fabian (1998) who find a tem- 
perature of 9.2910 60 keV (90%-confidence error bounds). 



For the X-ray surface-brightness fitting we used a 
ROSAT HRI image of A773 with an effective exposure of 
16518 s obtained on 13-15 April 1994 and analysed using 
standard ASTERIX routines. We calculate the ROSATHRI 
count rate, given our estimates of metallicity and Galactic 
column and with the K-correction appropriate to the red- 
shift of A773, to be 1.53(±0.08) x 10"^^ counts s"^ from a 
1 m^ cube of gas of electron density 1 m~* at the tempera^ 
ture of A773 and at a luminosity distance of 1 Mpc. 

We then fitted an ellipsoidal King profile to the X-ray 
image. Since the high spatial resolution of the HRI leads to 
a low count rate per pixel, we use Poisson rather than Gaus- 
sian statistics to fit for the measured count in each pixel. For 
Ci counts measured at position Xi, and for a mean rmmber 
f{xi\a) of counts predicted by the model given parameters a 
(such as core radius), the probability of obtaining a counts 
is 



and the most likely value of a can be obtained in a compu- 
tationally efficient way by maximizing 

InP(c|a) = ^^(ci In /(xi I a) — IncJ — f{xi\a)). 

i 

We fitted an ellipsoidal King profile to the HRI data 
with 6i and 02 as the perpendicular angular sizes in the 
plane of the image, assuming that the length along the 
line of sight is the geometric mean of the other two. We 
find 01 — 60" and 02 = 44", with the major axis at po- 
sition angle 16°, /3 = 0.64, and central electron density 
no = 6.80X 10^ /ijo^^V-^ where Ho = 50 hso km s"^ Mpc"\ 
Fig 5 shows the HRI image, the model, and the residual im- 
age with the best model subtracted. To assess the goodness 
of fit, we made 50 realisations of the image with the ap- 
propriate Poisson noise added, and calculated the mean and 
standard deviation of their Poisson likelihoods. The likeli- 
hood of the observed HRI image is 0.32 standard deviations 
from the mean; we therefore conclude that the fit is good 
and the cluster is well represented by a /3 model. 

There is a strong degeneracy in the fit between (3 and 
01,2; however this has little effect on the comparison with 
the SZ data and the derived value of Hq. Figure 6 shows the 
likelihood contours for the fit in the l3-0i plane, marginalised 
over no and using the best-fit value of the axial ratio (which 
is very well constrained) . Overlaid are the contours of pre- 
dicted mean observed SZ flux density on the shortest RT 
baseline. It can bo seen that despite the degeneracy between 
P and 01 , the range of SZ flux densities corresponding to the 
1-0" limits of the model fit is only ±3%. Since the SZ flux 
density varies as Hq^, this corresponds to a 6% error in Ho 
due to the model fitting. This lack of sensitivity to the (3-0 
degeneracy is characteristic of observations that are sensi- 
tive to spatial frequencies around the cluster core size (see 
eg Reese et al (2000)) and contrasts with the sensitivity to 
the model fitting of measurements that measure only lower 
spatial frequencies (eg Birkinshaw & Hughes(1994)). 

4 Ho ESTIMATION 

To measure Ho, we compared the real SZ data with a sim- 
ulation of the SZ effect from the X-ray gas model. We use 
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the expression of Challinor & Lascnby (1998) to provide 
a relativistic correction to the standard non-relativistic SZ 
expression; in the case of A773, the effect is to increase our 
estimate of the j/-parameter by 2.4%. We then simulated RT 
observations of the SZ effect due to the model gas distribu- 
tion and to compared these with the real source-subtracted 
RT visibilities on the same baselines, and adjusted Ho to 
get the best fit. Using our temperature of 8.7±0.7 keV we 
find Ho = ntu km s'^ Mpc'S assuming an Einstein-de- 
Sitter universe. The 1-a error quoted is that due solely to 
noise in the SZ data. For the best fit (3, 61,2 model, the 
corresponding central density no is 8.44 x 10^ m~^ and the 
central decrement 737 ± 85/iK. 

Graingc ct al (2002) consider at some length the con- 
tributions to error in the Hq determination from A1413. 
The situation in A773 is very similar. The dominant con- 
tributions to the error in Ho in A773 are ±16% from noise 
in the SZ measurement, ±12% from our estimation of the 
gas temperature and a likely error of ±14% from the un- 
certain hne-of-sight depth. This is obtained by considering 
the range of axial ratios of simulated clusters that is needed 
to reproduce the projected axial ratio distribution observed 
in clusters with redshift similat to that of A773 (Grainger 
2001). Clearly this estimate is rather uncertain for a sin- 
gle object, but can be significantly reduced by averaging a 
sample of clusters with random orientations. Table 2 shows 
the complete error budget, and the final l-cr error limits 
of Ho = 77ti5kms-iMpc-i if {Qu^Ma) = (1.0,0.0) and 
Ho = 85ti? kms"^ Mpc"^ if {n^, Oa) = (0.3, 0.7). 

5 CONCLUSIONS 

Using ASCA, ROSAT HRI, and RT observations of A773, 
we find: 

(i) there are eight radio sources detectable in the field of 
the cluster that we have removed from the data, which would 
otherwise contaminate the measurement of the SZ effect; 

(ii) the correlated fitting errors on the shape parameters 
P and 6 have negligable effect on the derived value of Ho, 
a feature characteristic of observations on the scale of the 
cluster core size; 

(iii) the estimated value of Ho is 77^\l kms~^ Mpc~^ if 
(!:!„„ Oa) = (1.0,0.0) or 85l 1? km s'^ Mpc"^ ff (nm,nA) = 
(0.3, 0.7) , where the 1-a error bars include estimates from 
the main sources of error — noise in the SZ data. X-ray tem- 
perature uncertainty, and uncertain line-of-sight depth. 
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Flux density //iJy RA offset /" dec offset /" 



163 
123 
55 
103 
228 
162 
238 
177 



-67 

-34 
-22 

-5 

35 

46 

90 
139 



-23 

9 

75 
-155 
-29 

96 

40 
-16 



Table 1 . Flux densities and positions of sources removed from the 
A773 field. Flux densities are apparent, ie not corrected for the RT 
primary beam response, and are all ±35 /iJy. Offsets are relative 
to the pointing centre of 09''17™51^91 + 51°43'32" (J2000). 
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Figure 1. Map of the longer basclinc-RT data (2kA upwards) 
(contours) plus a VLA 1.4 GHz observation (grcyscalc), jointly 
used to determine initial source positions for the source subtrac- 
tion procedure. The crosses show the positions of the subtracted 
sources listed in Table 1. The contour interval is 70 //Jy, dashed 
contours are negative; the greyscale runs from 0.3 (light) to 1.0 
mjy (dark). The RT data have not been corrected for the effect 
of the telescope primary beam (FWHM = 6'). 
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Figure 2. Map of the longer baseline data (2kA upwards) after 
source subtraction, showing no significant residual features. The 
crosses show the positions of the subtracted sources listed in Table 
1; the size of the cross is proportional to the flux subtracted. The 

grcyscalc runs from —200 (light) to 200 /ijy (dark). The data have 
not been corrected for the effect of the telescope primary beam. 
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Figure 3. X-ray and SZ images of A773. The greyscale is 
the ROSAT HRI image; contours axe the source-subtracted 
and CLEANed RT image. Contours axe in multiples of 
80 /iJy beam" ^ , dashed contours axe negative. The restoring beam 
(shown) is 152 X 119 arcsec, PA 4° and the data have not been 
corrected for the effect of the telescope primary beam. 
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Source of error 



Error contribution to Hq 



SZ measurement ±16% 

Line-of-sight depth ±14% 

Gas temperature ±12% 

Error in calculated X-ray emission constant ±6% 

Primary flux calibration ±2.5% 



Total (in quadrature) 



±25% 



Table 2. Error budget for Hq calculation; the dominant terms are the errors in the SZ measurement and the gas temperature. The 
errors are all multiplicative, ie the l-a limits on Ho are 77 x 1.25 = 96 and 77/1.25 = 62kms~^ Mpc~^. 





Figure 5. X-ray model fits, (left) ROSAT HRI image of A773. (centre) The best-fitting model X-ray surface brightness, (right) The 
difference between the data and the model, which is consistent with the Poisson noise in the data. The contour interval is 1 count per 
8" pixel. 
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Figure 4. The real part of the source-subtracted visibilities as a 
function of baseline in wavelengths (filled points with error bars), 
with the best-fitting SZ model (open points joined by lines). 
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Figure 6. Solid contours and grcyscalc: the likelihood of the fit- 
ted X-ray model as a function of /3 and 9i, marginalised over no, 
with the axial ratio and position angle fixed at their most likely 
values. The solid contours enclose 67% and 95% of the likelihood. 
Overlaid are dotted contours of the predicted mean flux on the 
shortest RT baseline for the given model parameters (using the 
most likely value of no); the contour interval is 40/iJy. Since the 
normalization of the SZ fiux relative to the X-ray is what con- 
trols the fit for h, this Figure shows the insensitivity of our Ho 
determination to the degeneracy in the X-ray model fit. 



© 0000 RAS, MNRAS 000, 000-000 



